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Abstract 

We compare theory and experiment in the Casimir force measurement between gold surfaces 
performed with the atomic force microscope. Both random and systematic experimental errors 
are found leading to a total absolute error equal to 8.5 pN at 95% confidence. In terms of the 
relative errors, experimental precision of 1.75% is obtained at the shortest separation of 62 nm 
at 95% confidence level (at 60% confidence the experimental precision of 1% is confirmed at the 
shortest separation) . An independent determination of the accuracy of the theoretical calculations 
of the Casimir force and its application to the experimental configuration is carefully made. Spe- 
cial attention is paid to the sample-dependent variations of the optical tabulated data due to the 
presence of grains, contribution of surface plasmons, and errors introduced by the use of the prox- 
imity force theorem. Nonmultiplicative and diffraction-type contributions to the surface roughness 
corrections are examined. The electric forces due to patch potentials resulting from the polycrys- 
talline nature of the gold films are estimated. The finite size and thermal effects are found to be 
negligible. The theoretical accuracy of about 1.69% and 1.1% are found at a separation 62 nm and 
200 nm, respectively. Within the limits of experimental and theoretical errors very good agreement 
between experiment and theory is confirmed characterized by the root mean square deviation of 
about 3.5 pN within all measurement range. The conclusion is made that the Casimir force is 
stable relative to variations of the sample-dependent optical and electric properties, which opens 
new opportunities to use the Casimir effect for diagnostic purposes. 

PACS numbers: 12.20.Fv, 12.20.Ds, 42.50.Lc, 05.70.-a 
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I. INTRODUCTION 



In the last few years the Casimir effect Q], which is a rare macroscopic manifestation 
of the boundary dependence of the quantum vacuum, has attracted much experimental 



). The spectrum of 



and theoretical attention (see monographs d |3|, |4( and reviews 
the electromagnetic zero-point oscillations depends on the presence of material bodies. In 
particular, the tangential component of the electric field vanishes on the surfaces of two 
parallel plates made of ideal metal (it is small if real metals are used) . This leads to changes 
in the zero-point oscillation spectrum compared to the case of free unbounded space and 
results in the attractive Casimir force acting normal to the surfaces of the plates. 



The Casimir effect finds many applications in quantum field t 
physics, elementary particle physics, gravitation and cosmology 
many measurements of the Casimir force have been performed 



reory, condensed matter 



15 




161 ] - Their results have already been applied in nanotechnology for the actuation of the 



novel microelectromechanical devices, based entirely on the modification of the properties 



of quantum vacuum 



171 ] , and for constraining predictions of extra-dimensional physics with 



low compactification scales 
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Most theoretical papers on the Casimir effect deal with idealized boundary conditions 
and perfectly shaped test bodies. Over the last 4 decades only a few have considered the 
corrections to the Casimir force such as due to the finite conductivity of the boundary metal 



h y, q, 



271 ] and nonzero temperature [2i 



distortions of the surface shape 
Comparison of the theory with the results of modern Casimir force measurements demands 
careful treatment of all these corrections. Both the individual corrections and their combined 
effect has to be evaluated (see Ref. (J for review). 

The quantification of errors and precision in the measurements and theoretical compu- 
tations of the Casimir force is crucial for using the Casimir effect as a new test for extra- 
dimensional physics and other extensions to the Standard Model. Nevertheless, there is no 
general agreement on the achieved levels of experimental precision and the extent of agree- 
ment between theory and experiment. In the literature a variety of measures to characterize 
the experimental precision is used and the extent of agreement between measurements and 
theory ranges from 1% Q, H H Q to 15% Q depending on the measurement scheme 
and configuration. Very often, the confidence levels and numerous background effects which 
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may contribute to the theoretical results are not considered. 

In the present paper we perform a reanalysis of the experimental data on the Casimir 
force measurements between Au surfaces jllj and make a comparison with theory. In doing 
so we carefully calculate the original experimental precision without relation to the theory, 
including the random absolute error at a 95% confidence level, and the absolute systematic 
error. The total absolute error of these Casimir force measurements in the experiment of 
Ref. is found to be equal to A tot « 8.6 pN at 95% confidence. This corresponds to 
approximately 1.75% precision at the closest separation a rs 62 nm (the 1% precision at 
the closest separation indicated in Ref ll| is obtained at 60% confidence). As a second 
step, the accu racy of the theoretical computations of the Casimir force for the experimental 
configuration |ll[ is determined. Special attention is paid to the possible sample-dependent 
variations of the optical tabulated data due to the presence of grains, contribution of the 
surface plasmons, and errors introduced by the use of the proximity force theorem. The 
influence of the surface roughness is carefully investigated inclu ding the nonmultiplicative 
contributions and recently discussed diffraction-type effects jsfll . Isij . The contribution of 
electric forces due to patch potentials resulting from the polycrystalline nature of the Au film 



is calculated for the experimental configuration [11| at different separations. The finite size 
and thermal effects are also considered and found negligible in the experimental configuration 
of Ref. [ll]. The conclusion reached is that at the present state of our knowledge the 
accuracies of theoretical computations in application to the experimental configuration of 
Ref. are achievable on the level of 1.69% at a separation z = 62 nm and 1.1% at a 
separation z = 200 nm. 

The paper is organized as follows. In Sec. II the experimental precision of the Casimir 
force measurements at different confidence levels is determined. Sec. Ill is devoted to the 
computations of the Casimir force with account of finite conductivity and grain structure of 
the metal layers. The role of roughness including the nonmultiplicative and diffraction-type 
effects is studied in Sec. IV. In Sec. V both traditional and alternative thermal corrections are 
discussed. Also the possible role of the electric forces due to the patch potentials and finite 
size effects are estimated. Sec. VI contains the final numbers on theoretical accuracy and 
the comparison of theory with experiment in the Casimir force measurement between two 
gold surfaces by means of an atomic force microscope ^l|. In Sec. VII the final conclusions 
and some discussion are provided. 
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II. EXPERIMENTAL PRECISION IN THE CASIMIR FORCE MEASURE- 
MENTS BETWEEN TWO GOLD SURFACES 



In Ref. precision measurements of the Casimir force between gold coated bodies, a 
plane plate and a sphere, were performed using an atomic force microscope. The Casimir 
force was measured by averagin g 30 scans over a surface separation region between 62-350 nm 
with 2583 points each (see Ref. [11 1 for all the details of the measurement procedure). In the 
analysis below we neglect data from 3 scans due to excessive noise and use the data from 
the rest n = 27 scans to find the quantitative characteristics of the experimental precision 
in the Casimir force measurements at different confidence levels. 

We start with the random error and calculate the mean values of the measured force at 
different separations z% within the region from 62 nm to 350 nm 

n 

j^(*) = ££jrw (!) 

k=l 

An estimate for the variance of this mean is determined by 

= E [*n*) - ^)] J - (2) 

^ ' k=l 

Calculations using the measurement data {F^ xp (zi)} show that s F (zi) do not depend sensi- 
tively on Zi. The largest value sp = 2.8 pN is taken below as an estimate for the variance of 
the mean force within the whole measurement range. 

According to Student's test for the truth of a hypothesis Q, if the inequality 

\F"*(z) - F(z)\ 

> *« = h-\cc \9) 

s F 2 

is fulfilled, the hypothesis that F(z) is the true value of the Casimir force at a separation 
z must be rejected at a given confidence level a (this is a two-tailed test as the deviations 
F(z) from F exp (z) in two directions are possible). Equivalently, if the inequality 

P^l < tl . la (4) 

Sp 2 

is fulfilled, the hypothesis that F(z) is the true value of the Casimir force should be accepted 
at a confidence level 8 = 1 — a. 

Usually in the tab.es for Student's .-distribution (see, e.g., Refs. Q Q) the values of 
t p = t p (f) are presented, where p — 1 — a/2 = (1 + 8)/2, f — n — 1 is the number of 



degrees of freedom, and n is the number of measurements (n = 27 in our case). Choosing 
(3 = 0.95 (hypothesis is true at 95% confidence) we obtain p = 0.975 and find from tables 
t p {f) = 2.056 H- Then from Eq. (gj) it follows 

\F exp {z) - F{z)\ < A rand F exp ee spt p (f) » 5.8 pN, (5) 

where A rand F exp is the random absolute error of the Casimir force measurements. If we 
consider (3 = 0.6 (hypothesis is true at 60% confidence), then p = 0.8 and t p (f) = 0.856. 
In this case the random absolute error of the Casimir force measurements is A rand F exp = 
2.8 x 0.856 pN^ 2.4 pN. Note that if one would like to have t p (f) = 1 or t p (f) = 2 (i.e. 
deviations of the true force value on either side of the mean not greater than one or two sp) , 
the confidence levels of (3 — 0.66 or (3 = 0.94, respectively, should be chosen for the number 
of measurements n = 27. 

Now let us consider the systematic error. The main contributions to the systematic error 
in the experiment of Ref. (llj are given by the error in force calibration F exp ~ 1.7 pN, 
by the noise when the calibration voltage is applied to the cantilever /\^ yst F ex P ~ 0.55 pN, 
by the instrumental sensitivity /\ s 3 yst F ex P ~ 0.31 pN, and by the restrictions on computer 
resolution of data A^ 8 F exp w 0.12 pN. The maximal value of the systematic error is given 
by 

4 

A syst F ex P = /\ s v st F^p ~ 2.7 pN. (6) 
1=1 

Finally, the maximum total absolute error of the Casimir force measurements in the exper- 
iment of Ref. jll| is equal to 

^jpexp _ ^rand pexp ^syst pexp ~ g_5 pN (7) 

at 95% confidence (to be conservative, the errors are added linearly rather than quadrat- 
ically). At 60% confidence the total absolute error of the Casimir force measurements is 
^pexp ^ 5_ipN_ These absolute errors with their confidence levels are valid within the 
whole measurement range from 62 nm to 350 nm. From Eq. (J3J) it follows that the true value 
of the Casimir force belongs to the confidence interval 

F exp {z) - AF exp < F{z) < F exp {z) + AF exp (8) 

with a chosen confidence probability. 
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Another important characteristic of the experimental precision is the relative error of 
the Casimir force measurements SF exp (z) = AF exp /F exp (z) which is evidently separation- 
dependent. At the shortest separation z = 62 nm the value of the mean force is F exp = 
485.8 pN which leads to a relative error of 5F exp (z) « 1.75% computed at 95% confidence. 
If we restrict ourselves with a 60% confidence, the relative error of the Casimir force mea- 
surements at the shortest separation 5F exp (z) = 5.1/485.8 ~ 1% is obtained as was indicated 



in Ref. 



ll\ without the detailed analysis of the confidence levels. If we choose 95% confi- 



dence, the relative errors of the Casimir force measurements at separations 70 nm, 100 nm, 
and 200 nm are, respectively, 2.46%, 5.9%, and 37.3%. At 60% confidence the relative errors 
of the Casimir force measurements at the same separations are 1.47%, 3.5%, and 22.4%, 
respectively. 



III. CALCULATION OF THE CASIMIR FORCE INCLUDING THE FINITE 
CONDUCTIVITY AND GRAIN STRUCTURE OF GOLD LAYERS 

For the configuration of a large sphere of a radius R above a plate the Casimir force can 
be obtained by means of the Lifshitz formula, derived originally for two parallel plates [34] , 
along with use of the proximity force theorem 

f- D roo POO 

F c {z) = / fc A / di {In [1 - r}(£, k^e^] 
^ Jo Jo 

+ ln[l-ri(£,£; ± )e- 2 ^}. (9) 
Here the reflection coefficients for two independent polarizations are given by 



q — k 
q + k 



q 2 = k'i + ^, k * = kl + e(i0^ (l0 ) 
c c 



e{uS) is the dielectric permittivity of the gold layers on the sphere and the plate, and z is the 
closest separation distance between them. The thickness of gold coatings, used in Ref. [ill 
(86.6 nm), is much greater than the skin depth of the electromagnetic oscillations for all 
frequencies which make a significant contribution to the computation of the Casimir force. 
This allows one to use the properties of the bulk gold in all computations of the Casimir 
force. 
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The accuracy of Eq. Q is restricted by the accuracy of the proximity force theorem, which 
is, however, very high for the experimental parameters of Ref. The error, introduced by 



371 ] . Taking into account the large value 



the proximity force theorem, is less than z/ R l36l 
of sphere radius R = 95.65 /im, used in Ref. the upper limit of this error is 0.06% at the 
shortest separation z = 62 nm and 0.2% at separation z = 200 nm (note that in Ref. [36j the 
Casimir force for the configuration of a sphere above a plate was precisely computed on the 
basis of the first physical principles which makes it quite reliable as a test of the proximity 
force theorem) . 

In Refs. the computations of the Casimir force were performed using Eqs. (jUJ), ffTTTft 

and optical tabulated data for gold j3] (note that the transition coefficient from energies 
to frequencies is given by leV= 1.52 x 10 15 rad/s). The imaginary part of the dielectric 
permittivity, obtained using the complex refractive index from the Tables j^j], was used 
to compute the dielectric permittivity along the imaginary frequency axis by means of the 
dispersion relation. At uj < 1.9 x 10 14 rad/s, where the tabulated data are not immediately 
available, they were usually obtained (see, e.g., by the extension from the region of 

higher frequencies by means of the Drude dielectric function 

e{u) = 1 - , ^ , , , (11) 
uj(uj + 27) 

where the plasma frequency for Au is uo p = 1.37 x 10 16 rad/s, and 7 = 5.32 x 10 13 rad/s is 
the relaxation parameter describing the non-elastic electron-phonon interaction (note that 
in the frequency region under consideration 7 <C uj). This procedure was used to calculate 
the Casimir force including the effect of finite conductivity corrections of gold (see a few 
examples of the calculations in Sec. IV, Table II and comparison between experiment and 
theory in Sec. VI). Later in this section we discuss the influence of possible sample to 
sample variations of the optical tabulated data on the values of the Casimir force and the 
applicability region of Eqs. ©, (|TU|) involving the dielectric permittivity depending only on 
frequency. 

First, we would like to note that in the separation region 200 nm < z < 350 nm the 
computational results obtained by Eq. Q combined with the optical tabulated data, are 
almost exactly those obtained by the substitution into Eq. (0) of the plasma dielectric 
function for the metal 

e{uj) = l-f 2 . (12) 
8 



In fact, both computations lead to results differing by less than 0.5% within the mentioned 
separation region. What this means is that the real part of e depending on only u p determines 



the total value of the Casimir force in this region. The value of u p = 2y/nNe/ y/rnF, where 
TV is the density of conduction electrons, m* is their effective mass, is determined by the 
preperties of the elementary cell. It cannot be influenced by properties of sample such as 
the crystallite grain size or the presence of a small concentration of impurities. This is the 
reason why the sample to sample variations of the optical tabulated data cannot influence 
the value of the Casimir force at separations z > 200 nm. 

In the separation region 62 nm < z < 200 nm there are significant deviations depending on 
whether the Casimir force (jHJ) is calculated using the optical tabulated data or by use of the 
plasma dielectric function ()12|1 . In fact, in this separation region the small imaginary part of e 
is influential and should be taken into account. There is enough tabulated data in the optical 
Tables to compute the Casimir force, so that it is not necessary to use any extension of data. 
Note that the characteristic frequency corresponding to the largest separation z = 200 nm 
is u c = c/(2z) = 7.5 x 10 14 rad/s (i.e. tabulated data for frequencies several times smaller 
are available). At the same time, the characteristic frequency corresponding to the shortest 
separation is u c = 2.42 x 10 15 rad/s <C u p , so that the region under consideration belongs 



to that of infrared optics [4l|. Within the region (62-200)nm one may expect some small 
dependence of the optical tabulated data on the size of the grain, presence of impurities 
etc. If this is indeed the case, the use of the tabulated data, which are not relevant to 
the particular samples used in experiment, might lead to the errors in computation of the 
Casimir force (|9*|). 

To investigate this possibility, we consider the pure imaginary part of the dielectric per- 
mittivity in the region of the infrared optics given by [42] 

2 

Ime = (13) 

where v is the relaxation parameter at high frequencies in the region of infrared optics (note 
that it does not coincide with the relaxation parameter 7 of the Drude model which 
describes the volume relaxation in the region of the normal skin effect). The value of v 
is determined by the processes of elastic scattering of the electrons on impurities, on the 
boundary surfaces of the metal and of the individual grains, and on other electrons Q] • 
The scattering of electrons on phonons also contributes to the value of v. However, the 



frequency of the electromagnetic field is so high that hu 3> &bTd, where To is the Debye 
temperature, so the frequency of the electron-phonon collisions is the same as it is at T = To 



42| . It is important to note that of all the above processes, only the contribution of the 
electron-electron collisions to v is frequency-dependent (and increases as uj 2 ). 

The main sample to sample dependence of the parameter v is determined by the sizes 
of grains and the density of impurities. To calculate this dependence we use the following 



formula for the relaxation parameter in the region of infrared optics 



42. 



44] 



^) = l + ^-^(c 1 -c 2 i,), (15) 



v = uj p ^ci + c 2 ^j ■ (14) 

This formula leads to an approximate representation of the dielectric permittivity of Au 
along the imaginary frequency axis given by 

i «i( e 

where c\ = 0.0039, c 2 = 1.5. It is easily seen, that the substitution of Eq. (|T5j) into Eq. 
leads approximately to the same result as the use of the optical tabulated data. The errors 
due to use of Eq. (fTHjl in Eq. instead of the optical tabulated data at separations 62 nm, 
70 nm, 100 nm, and 150 nm, are 0.45%, 0.23%, 0.09%, and 0.04%, respectively. 

Eq. ()15j) gives the possibility to estimate the influence of the sizes of grains in the poly- 
crystalline metal film in the experiment of Ref . II 1 1 on the value of the Casimir force Q. 
For this purpose, the experimental data of Ref. |45j are used where the reflectance R of Au 
films is measured as a function of the characteristic sizes of the grains. 

The analysis of the atomic force microscopy images (like the one in Fig. 1 but on 1 x 1 /xm 2 
area) shows that the mean size of grains in Ref. [ll| is about 90 nm (the sizes of the typical 
grains are 77 nm, 103 nm, 94 nm, 68 nm, 88 nm, 121 nm etc). According to Ref. 45j, the 
largest deviations of the reflectance from the one given by the tabulated data 40j, takes 
place at shorter wavelengths. The shortest separation of z = 62 nm in the experiment ^j] 
corresponds to the characteristic wavelength A c = 2ttc/uj c = 4irz pa 780 nm. For the films 
containing grains of 45 nm size (the largest ones studied in Ref. j^J) the reflectance at 
A ~ (750 — 800)nm is 0.8% less than the one calculated from the tabulated data (notice 
that for smaller grains the difference of the reflectance obtained from the tabulated data 
is greater). Taking into account that the reflectance in the region of the infrared optics is 
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given by 



i?=l-4Re^ = — , (16) 

we find that the new value for the coefficient c\ in Eqs. (114)1 . ()15j) due to grains of 45 nm size 
is c\ = 0.0059. Substituting the approximate Eq. (|TH)) (with c\ instead of C\) into Eq. Q, 
one finds the values of the correction factor to the Casimir force fj^ = 0.439 at z = 62 nm 
and fjf = 0.465 at z = 70 nm. Comparing this with the results of the same approximate 
computations using c\ (rj^ = 0.441, respectively, r/^ = 0.467), one can conclude that the 
grains of 45 nm size lead to less than 0.5% decrease of the Casimir force magnitude. Note 
that this is in fact the upper bound for the influence of crystallite grain size on the Casimir 
force in the experiment of Ref . ^l] , as the actual sizes of grains in were two times greater 
than 45 nm. 

The above calculations of the Casimir force including the effect of the real properties of Au 
films were performed on the basis of the Lifshitz formula (jUJ), which does not take into account 
the effects of spatial nonlocality (wavevector dependence of the dielectric permittivity). 
These effects may influence the Casimir force value in the region of the anomalous skin 
effect which is important for large separations z > 2.36 //m |4jJ], a region not relevant to 
the experiment of Ref. ll| . Another separation region, where nonlocality may lead to 
important contributions to the van der Waals force, is z < X p /(Air) « 10.9 nm (A p is the 
plasma wavelength) which corresponds to u c > u p |48j . Such high characteristic frequencies 
lead to the propagation of surface plasmons. The effect of the surface plasmons, however, 
does not contribute in the experiment of Ref. as the largest characteristic frequency 
there, calculated at z — 62 nm, is 5.7 times less than uo p [notice that the frequency region 
(5u c , I0u> c ) contributes only 0.19% of the Casimir force value at separation z = 62 nm]. The 
contribution of the surface plasmon for Au of about 2% at a separation z = X„ = 137 nm, 



obtained in recent Ref. 49], is explained by the use in [49] of the spatially nonlocal dielectric 
permittivity in the frequency region of infrared optics where it is in fact local 
(we would like to point out that at the separation z = X p the characteristic frequency of the 
Casimir effect is equal not to uj p , as one might expect, but u p / (47r)]. As a result, the surface 
plasmons do not give any contribution to the Casimir force in the experimental configuration 



of Ref. 



ui- 
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IV. SURFACE ROUGHNESS CORRECTION TO THE CASIMIR FORCE AND 
ITS CALCULATION USING DIFFERENT APPROACHES 



It is well known that surface roughness corrections may play an important role in Casimir 



force calculations at separations less than 1 /zm (|. At the shortest separations, the rough- 
ness correction contributes 20% of the measured force in experiments of Refs. 8|, [l5j, llq . 
In the experiment of Ref. [ill ], however, the roughness amplitude was decreased and the 
roughness contribution was made less than 1% of the measured force even at shortest sep- 
arations. To obtain this conclusion the simple stochastic model for the surface roughness 
and the multiplicative approach to take into account different corrections were used. Here 
we obtain more exact results for the contribution of surface roughness to the Casimir force 
taking into account both nonmultiplicative and correlation effects. 

The topography of the Au coatings on the plate and sphere was investigated using an 
atomic force microscope. A typical 3-dimensional image resulting from the surface scan 
of 15 /im x 15 /im area is shown in Fig. 1. As seen in this figure, the roughness is mostly 
represented by the stochastically distributed distortions with the typical heights of about 
2-4 nm, and rare pointlike peaks with the heights up to 16 nm. In Table I the fractions Vi 
of the surface area, shown in Fig. 1, with heights hi are presented (i = 1, 2, . . . , 17). These 
data allow one to determine the zero roughness level Hq relative to which the mean value of 
the function, describing roughness, is zero (note that separations between different bodies 
in the Casimir force measurements are usually measured between the zero roughness levels 
3): 

17 

Y,(H -h i )v i = 0. (17) 

i=l 

Solving Eq. (fTTj). one obtains H rs 2.734 nm. If the roughness is described by the regular 
(nonstochastic) functions Af(x,y), where \f(x,y)\ < 1, for the roughness amplitude it 
follows A = h™ ax -H = 13.266 nm. 

In the framework of the additive approach the values of the Casimir force including the 
effect of finite conductivity F c (z), obtained in Sec. Ill, Eq. ©, can be used to calculate the 
effect of roughness. For this purpose, the values of F c should be geometrically averaged over 
all different possible separations between the rough surfaces weighted with the probability 
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of each separation 



HQ 
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Fc,t{z) = ViVjF c (z + 2H -hi - hj 
m=i 



(18) 



Note that Eq. (jTSj) is not reduced to a simple multiplication of the correction factors due to 
finite conductivity and surface roughness but takes into account their combined (nonmulti- 
plicative) effect. 

An alternative method of calculating the corrections due to the stochastic surface rough- 
ness was used in Ref. According to the results of Ref. 5(|, the Casimir force between 
a plate and a sphere made of ideal metal and covered by a stochastic roughness with an 
amplitude A st is given by 



F r (z) = F (z) 



1 + 6 



A 



si 



+ 45 



.4 



si 



(19) 



where F (z) = —ii 3 hcR/(360z 3 ) is the Casimir force between perfectly shaped plate and 
sphere of radius R. Then the Casimir force including both the finite conductivity of the 
boundary metal and surface roughness can be calculated as 



Fc{z) 



01 — 

z 



z 



(20) 



i.e. by means of the multiplicative procedure. 

The variance of the random process describing the stochastic roughness is found by the 
formula 

17 

fc^iHo-hifvi. (21) 



i=l 



Using the data from Table I, one obtains the values for variance 5 s t ~ 0.837 nm and for the 
amplitude of a random process A st = V~25 s t ~ 1.18 nm. This value is slightly larger than 
the one obtained in Ref. 1]] on the basis of less complete data on roughness topography. 

Now we are in a position to compare the contribution of the surface roughness computed 
by Eq. ([THjl . taking into account the combined effect of the roughness and finite conductivity, 
and by the multiplicative procedure of Eq. (|20|). In Table II the results for the correction 
factors r] c = F c f F , r\ T = F r /F , r\ c ^ r = F c ^ r /F , and r}™ T = r] c r] r are presented at the shortest 
separations z = 62 nm, 70 nm, 80 nm, and 90 nm, where the roughness corrections play some 
role. As is seen from Table II, both approaches lead to practically coincident results for the 
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roughness correction factors due to the combined effect of finite conductivity and surface 
roughness. This means that for such small roughness as in Ref. [11] the multiplicative pro- 



cedure is quite satisfactory (for 
contributions may be essential 



arger roughness amplitudes, however, the nonmultiplicative 
H lid]). Note also that for A st w 1.18 nm the fourth order 
term in Eq. (|2(J|) practically does not contribute even at shortest separations and can be 
neglected as was done in Ref. 

Both Eqs. (fTSj) and pUj) used above are based on the approximation of additive summation 
and do not take into account the diffraction-type effects which arise in the case of roughness 
described by the periodic functions with small periods A < z (^J or by the stochastic 
functions with small correlation length 3l( . To estimate the value of the correlation length 
in our case, we consider the set of cross sections of the roughness image shown in Fig. 1. 

In Fig. 2 two typical cross sections are presented, one at fixed x (a) and the other one 
at fixed y (b). We have performed the Fourier analysis of the functions, as in Figs. 2,a,b, 
along the lines of Ref. j22|- It was found that the Fourier harmonics, giving the major 
contribution to the result, are characterized by significantly greater periods than the mean 
distance between the neighbouring peaks in Figs. 2,a,b which is equal, approximately, to 
180 nm. 

To obtain an estimate for the upper limit of the contribution of the diffraction-type effects 
in the above roughness analysis, we use the correlation length l corr = 200 nm (slightly larger 
than the mean distance between peaks) and consider the periodic function with this period 
(clearly, the diffraction-type effects are greater for a periodic function with a period l corr 
than for the random function with a correlation length lcorr)- With this the diffraction-type 
effects can be computed in the framework of the functional approach developed in Ref. [30] . 
At a shortest separation z = 62 nm one obtains z/l corr ~ 0.31. Then for the coefficient c corr 
in the expression 

A st x ' 2 



r/r r = l + 6cw — , (22) 

n 

taking the diffraction-type effects into account, from Fig. 2 of Ref. |3J| it follows c corr ~ 1.1. 
As a result, using the upper limit for the contribution of the diffraction effects one obtains 
, rf° rr pa 1.0024, i.e. only 0.02% difference with the value of rj r in Table II obtained by 
neglecting the diffraction effects. At larger separations the diffraction effects lead to larger 
contribution to the roughness corrections. For example, at a separation z = 90 nm we 
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have z/lcorr ~ 0.45, c corr m 1.28, and rf° rr m 1.0013, i.e. 0.03% difference with the result 
of Table II. At larger separations, however, the roughness correction itself is even more 
negligible than at the shortest separations. 

To conclude, the surface roughness contribution in the experiment of Ref. |h| does not 
exceed 0.24% of the Casimir force at the shortest separation z = 62 nm. The diffraction-type 
effects, which were not taken into account in Eqs. ()18j) . (|20p . are shown to contribute less 
than one tenth of this result. 



V. CONTRIBUTIONS OF THE THERMAL CORRECTIONS, RESIDUAL ELEC- 
TRIC FORCES AND FINITE SIZES OF THE PLATE 

Although the experiment of Ref. ^j] was performed at room temperature T = 300 K, all 
the above computations were done at zero temperature. The thermal Casimir force F c (z, T) 
is given by Eq. Q where integration in continuous £ is changed to a summation over the 
discrete Matsubara frequencies £j = 2nkBTl/h according to 



CO 



AC 2nk B T ^ / 



h 

1=0 

leading to the Lifshitz formula for the thermal Casimir force. Here prime refers to the 
addition of a multiple 1/2 near the term with I = 0. When T — ► 0, F c (z,T) — > F c (z, 0) = 
F c (z), where F c (z) is given by Eq. (J0J). 

The magnitude of the relative thermal correction to the Casimir force can be computed 
by the formula 

F c (z,T) - F e (z) _ A T F c (z,T) 
5tF ^ Z ' T) = JM = FJz) ■ (23) 



Recently, there has been extensive discussion in 
procedure for the thermal Casimir force F c (z,T) 



iterature on the correct calculation 
In Refs. |5jJ, k52| the dielectric 



HQ 

permittivity of the plasma model ()12|) was substituted into the Lifshitz formula for F c (z, T). 
This approach, which was later called "traditional" 0], leads to the thermal corrections 
A l rfF c , 5 l r£F c . It is consistent with thermodynamics and agrees with the limiting case of the 
ideal metal. In the region of infrared optics the same results were obtained in the framework 
of the impedance approach which does not consider the fluctuating electromagnetic field 
inside the metal and takes into account the realistic properties of the metal by means of the 

47, 53 1 . Within the separation distances of Ref. Illl, the 



Leontovich boundary condition 



53(. Within the separation distances of Ref. jll|. 
15 



traditional thermal corrections are very small. As an example, at a separation z = 100 nm 
and T = 300 K one has 6!fF c w 0.007%, and 5%F C pa 0.03%, 0.1% at separations z = 200 nm, 
respectively, 300 nm js^J (to compare, in the case of ideal metals the same corrections, found 
in the framework of the thermal quantum field theory, are equal to 0.003%, 0.024%, and 
0.08%, respectively, i.e. the results for real metals approach the results for ideal ones with 
the increase of separation 0, El)- Thus, the traditional thermal corrections are negligible 



in the measurement range of experiment 



111 ] (the contribution of the relaxation processes 



to the magnitude of these corrections, which can be computed by taking into account the 
small real part of the surface impedance, is much less than the corrections). 

Alternatively, in Refs. 55^ |56[ the dielectric permittivity of the Drude model (jll|) was 
used to calculate F c (z,T). In this approach, there is no continuous transition between the 
cases of real and ideal metal. At the high temperature limit, the Casimir force between real 
metals was found equal to one half of the result obtained for the ideal metal (independently 
of how high the conductivity of real metal is). The thermal corrections, computed in the 



tal is). r 

0, Q, 



framework of the alternative approach 3], are quite different from those obtained 

from the traditional approach. To find the magnitude of these corrections [54], one should 
substitude into Eq. (|2*3l) 

A T F c (z,T) = A?F c (z,T) « A*F c (z,T) - / ydyki [l - r 2 ± (0, y) e -»] , (24) 

° a Jo 

where r^(0, y) is obtained by the substitution of Eq. ()12|) into Eq. (|10p. After calculations, 
one obtains that the alternative relative thermal correction increases from 5^F r , w 1.1% 



and 1.3% at separations z = 62 nm, respectively, 70 nm to 5^ F c pa 8% at a separation 
z = 350 nm. 

Another alternative thermal correction suggested in literature is also based on the 
substitution of the Drude dielectric function ()11|) into the Lifshitz formula for F r (z, T) but 
with a modified zero-frequency contribution for the perpendicular mode (in Ref. [57j this 
contribution is postulated to be of the same value as for an ideal metal). The alternative 

n n 

thermal correction of Ref. [5j| is given by |54j 

A T F c (z,T) = A { ?F c { Zl T) ^ A?F c (z,T) + -j-^ / ydyln (l - c ~v) . (25) 

(2) 

As a result, the relative alternative thermal correction of this kind takes values 5 T F c pa 
(2.1 — 2.2)% at all separations from z = 62 nm to z — 350 nm, i.e. slightly larger than the 
experimental precision at the shortest separations. 
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55 



561 ] and of 



As was shown in Ref. j^, both alternative thermal corrections of Refs. 
Ref. (5?! are not consistent with thermodynamics leading to the violation of the Nernst heat 
theorem. Recently they were found to be in disagreement with the precision measurement 
of the Casimir force using a microelectromechanical torsional oscillator [15]. In Sec. VI we 
will discuss the influence of the alternative thermal corrections on the comparison of theory 
and experiment in the Casimir force measurement of Ref. • 

In the rest of this section we discuss the probable contribution of the residual electric 
forces and the finite sizes of the plate on the Casimir force. As was noted in Ref. [llj], the 
electrostatic force due to the residual potential difference between the plate and the sphere 
has been lowered to negligible levels of -C 1% of the Casimir force at the closest separations. 
In recent Ref. (5^ it was argued, however, that the spatial variations of the surface potentials 
due to the grains of polycrystalline metal (the so called " pat ch potentials" ) may mimic the 
Casimir force. Here we apply the general results of Ref. 59] to the experiment of Ref. Q] 
and demonstrate that the patch effect does not make significant contributions. 

According to Ref. j^], for the configuration of a sphere above a plate the electric force 
due to random variations in patch potentials is given by 

_ 4ne o* v R /*« k^dk 

r v\ Z ) ~ k 2 _t,2 / „; nh t, 7 ' 

where a v is the variance of the potential distribution, fc max (fc mm ) are the magnitudes of 
the extremal wavevectors corresponding to minimal (maximal) sizes of grains, and Eq is 
the dielectric permittivity of free space. The work functions of gold are V\ = 5.47 eV, 
Vi = 5.37eV, and V3 = 5.31 eV for different crystallographic surface orientations (100), 
(110), and (111), respectively. Assuming equal areas of these crystallographic planes one 
obtains 

1 



i=i 



1/2 



80.8 mV. (27) 



Using the atomic force microscopy images discussed in Sec. Ill, the extremal sizes of 
grains in gold layers covering the test bodies were determined A min ~ 68 nm, and A max ~ 
121 nm. This leads to /c max ~ 0.092 nm -1 and fc mm pa 0.052 nm -1 . Note that these grain 
sizes are of the same order as the thickness of the film. The computations by Eq. (|2*6*|) 
using the above data lead to the "patch effect" electric forces F p /R pa —1.15 x 10~ 8 N/m 
and —1.25 x 10~ 10 N/m at separations z = 62 nm and z = 100 nm, respectively. Comparing 
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the obtained results with the values of the Casimir force at the same separations {F c /R pa 
—5.06 x 10~ 6 N/m, respectively, —1.48 x 10 _6 N/m), we conclude that the electric force 
due to the patch potentials contributes only 0.23% and 0.008% of the Casimir force at 
separations z = 62 nm, respectively, z = 100 nm (at a separation z = 200 nm the patch effect 
contributes only 7 x 10 _7 % of the Casimir force. So a rapid decrease of the contribution of 
the electric force with an increase of a separation is explained by the exponential decrease of 
the integral in Eq. (j2T?j) if one substitutes the physical values of the integration limits based 
on the properties of gold films used in experiment of Ref. 11]. 

A more important contribution of the patch electric forces may be expected in the scan- 
ning of a sharp tip of the atomic force microscope at a height of about 10-20 nm above a 
plate. In this case the electric forces are comparable with the van der Waals forces compli- 



cating the theoretical interpretation of force-distance relations [6 

Let us finally estimate the theoretical error caused by finiteness of the plate used in the 
experiment of Ref. jll| . Eq. (pj) was derived for the plate of infinite radius. In fact, the 
radius of the plate used in the experiment is L = 5 x 10~ 3 m. In this case, the Casimir 



force can be obtained by the formula 



F?»{z) = F c (z)P(z) 



i-ifi 



F c (z). (28) 



R 3 \ y/l + L 2 /R 2 

To calculate /3(z), we put z = 350 nm (to make this factor maximally distinct from unity) 
and obtain 

z 3 R 3 

(3(z) ^ 1 - 8 — — pa 1 - 2.2 x 10~ 17 , 

i.e. the finiteness of the plate size is too small to give any meaningful contributions to the 
Casimir force. 



VI. THEORETICAL ACCURACY AND COMPARISON OF THEORY AND EX- 
PERIMENT 

Now we are in a position to list all the sources of errors in the theoretical computation 
of the Casimir force F cr given by Eq. ()18j) . to find the final theoretical accuracy and to 
consider the comparison of theory and experiment. 

The main error, which arises from Eqs. 0, ()18jl when one substitutes the experimental 
data, is due to the errors in determination of separation z and sphere radius R. In terms of 
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dimensionless variables, F CjT .(z) is proportional to R and inversely proportional to z 3 which 
results in 

^ = ^k^ +3 ^. (29) 

In Ref. the absolute separations were determined by means of the electric measurements 
which allowed the determination of the average separation distance on contact zq with the 
absolute error A el z w 1 nm. Contrary to the opinion expressed in Ref. [62], this error 
should not be transferred to all separations leading to rather large contribution to 5^F Cjr 
of about 3A el zo/z « 4.8% at the shortest separations. Note that we are comparing not one 
experimental point with one theoretical value, but the experimental force-distance relation 
with the theoretical one computed on the basis of a fundamental theory. Thus, an additional 
fit should be made, with z$ as a fitting parameter changing within the limits (zq — A el zo, zq + 
A el zo) to achieve the smallest root mean square (r.m.s.) deviations between experiment and 
theory 

f M A V2 

aM = | m S ^ Fc ' r{Zi) ~ ) ' (30) 

where F exp (zi) was defined in Eq. ({I]). F c>r (zi) were computed by Eq. (j!8|) . and M is the 
number of experimental points under consideration. If, as usual, we consider two hypotheses 
as equivalent when they lead to the r.m.s. deviations differing for less than 10%, this results 
in decrease of the error in determination of absolute separations up to Az ~ 0.15nm. 

It is important to underline that the verification of the hypothesis is performed within 
different separation intervals (i.e. the total number M = N = 2583 experimental points 
within the whole separation range from 62 nm to 350 nm, M = 1270 points belonging to 
the interval 62-210 nm, and M = 600 points at separations less than the plasma wavelength 
X p = 136 nm). The above value for Az is almost one and the same in all the separation 
intervals. The obtained values of the r.m.s. deviations between theory and experiment are 
o"at ~ 3.4 pN, <7i27o ~ 3.2 pN, and o"6oo ~ 3.8 pN. These values are rather homogeneous 
demonstrating good agreement between theory and experiment independently of the chosen 
separation region. 

The radius of the sphere was measured more precisely than in Ref. [l]J with a result 
2R = 191.3 ± 0.3 fim. Using this together with Az = 0.15 nm, one obtains from Eq. ()29)1 at 
the shortest separations 5^F C:r m 0.88%. 

Now let us list the other contributions to the theoretical error of the Casimir force com- 
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putations at the shortest separation z = 62 nm and indicate their magnitude. According to 
the results of Sec. Ill, the sample to sample variations of the optical tabulated data may lead 
to the decrease of the Casimir force magnitude for no more than 8^ 2 'F Cjr fa 0.5%. The use of 
the proximity force theorem at z = 62 nm leads to very small error of about S^F cr fa 0.06% 
(Sec. III). The corrections due to the surface roughness are already incorporated in the theo- 
retical expression (|18|) but the diffraction-type effects may contribute up to 8^'F Cir fa 0.02% 
(Sec. IV). The effect of electric forces due to the patch potentials contribute a maximum 
8^F C:r fa 0.23% at the shortest separations, as was shown in Sec. V. The corrections due to 
the surface plasmons and finite size of the p 
and experimental configuration used in Ref. 



ate are negligible for the separation distances 
111 (see Sees. Ill and V). 



Special attention should be paid to the thermal corrections to the Casimir force. Accord- 
ing to the results of Sec. V, the contribution of the traditional thermal correction at the 
shortest separation is negligible. At larger separations it may be incorporated into the theo- 
retical expression for the force. As to the alternative thermal corrections of Refs. 
which contribute of about (1-2)% of the Casimir force at the separation z = 62 nm, they 
have been already ruled out both experimentally and theoretically (see Sec. V). If we would 
include any of these corrections into the theoretical expression for the Casimir force, this 
results in the increase of the r.m.s. deviation between theory and experiment which cannot 
be compensated by shifts of the separation distance in the limit of error A el z . In view of 
the above, we exclude the contributions from these hypothetical corrections from our error 
analysis. 

The upper limit for the total theoretical error at a separation z = 62 nm can be found by 
the summation of the above contributions 

5 

8F c , r = 6(l)F c,r « 1.69%, (31) 

i=l 

which is a bit more accurate than the total experimental relative error at the shortest 
separation equal to 1.75% at 95% confidence (see Sec. II). 

Note that with the increase of separation the experimental relative error quickly in- 
creases to 37.3% at a separation z = 200 nm. At the same time, the theoretical error is 
slowly decreasing with increasing separation. Thus, at z — 200 nm the above contributions 
to the theoretical error of the Casimir force computations take values 8^F Cj1 . fa 0.38%, 
5^F C:r fa 0.5%, 5^F C:T fa 0.21%, 5 (4) F Cjr fa 0.026%, 5&F c>r fa 0.000%. As a result, the total 
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theoretical error at z — 200 nm is 5F c ^ r pa 1.1%. 

The obtained results demonstrate very good agreement between theory and experiment 
within the limits of both experimental and theoretical errors. 

VII. CONCLUSIONS AND DISCUSSION 

In this paper we have performed a detailed comparison of experiment and theory in the 
Casimir force measurement between the gold coated plate and sphere by means of an atomic 
force microscope [llj|. The random error of the experimental values of the Casimir force was 
found to be A rand F exp pa 5.8 pN at 95% confidence (at 60% confidence the value A rand F exp pa 
2.4 pN was obtained). Together with the systematic error A syst F exp pa 2.7pN, this leads to 
the total absolute error of the Casimir force measurements in Ret. Q AF"' ss 8.5 pN at 
95% confidence. In terms of the relative errors, the experimental precision at the shortest 
reported separation is equal to 1.75% (1%) at 95% (60%) confidence level. 

In order to find the theoretical accuracy of the Casimir force calculations in the experi- 



mental configuration of Ref. many corrections to the ideal Casimir force were analysed. 
The correction due to the finite conductivity of gold was computed by the use of the opti- 
cal tabulated data of the complex refractive index. The results were compared with those 
computed by the use of the plasma dielectric function and found to coincide for the surface 
separation range (200-350)nm. At shorter separations the use of the optical tabulated data 
more accurately represents the dielectric properties. A special model was presented, which 
allows to take into account the sample to sample variations of the optical tabulated data 
due to the sizes of grains and impurities. It was shown that the error introduced by the 



grains of 45 nm size (even smaller than those in the experiment of Ref. 11 1) does not exceed 
0.5% of the Casimir force. The influence of the surface plasmon in the separation region of 
the experiment |ll( was found to be negligible. 

The surface roughness of the test bodies, used to measure the Casimir force, was carefully 
investigated by means of the atomic force microscope with a sharp tipped cantilever instead 
of a large sphere. The obtained profiles of roughness topography allowed calculation of 
the roughness corrections to the Casimir force in the framework of both multiplicative and 
nonmultiplicative approaches. The minor differences in the size of the effect are found only 
at the shortest separation. The correlation length of the surface roughness on the test bodies 
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was estimated and the diffraction-type effects were computed. At the shortest separation the 
roughness correction contributes 0.24% of the Casimir force with account of diffraction-type 
effects (and 0.22% with no account of diffraction). 

The electric forces caused by the spatial variations of the surface potentials due to the size 
of grains were investigated for the experimental configuration of Ref. [llj. They were shown 
to contribute 0.23% of the Casimir force at the shortest separation, and this contribution 
quickly decreases with an increase of separation. Several other effects (such as thermal 
corrections, corrections due to the finiteness of the plate and due to the deviation from the 
proximity force theorem) were investigated and found to make only negligible contributions. 

The final theoretical accuracy of the Casimir force calculations in the experimental con- 



figuration of Ref. 11 1 is 1.69% at the shortest separation z = 62 nm and 1.1% at a separation 
z = 200 nm. In the limits of both experimental and theoretical errors, very good agreement 
between theory and experiment was demonstrated characterized by the r.m.s. deviation of 
about 3.5 pN (less than 1% of the measured force at a shortest separation) which is almost 
independent of the separation region and the number of the experimental points. The above 
analysis does not support the conclusion of Ref. J3] that to achieve the experiments on the 
Casimir effect with 1% precision it is necessary to measure the separation on contact Zq with 
atomic precision. 

The obtained results demonstrate that in fact the Casimir force is more stable, than one 
might expect, to some delicate properties of the metallized test bodies like the variations 
of the optical data, patch potentials, correlation effects of roughness etc. These properties 
may change from sample to sample leaving the basic character, and even the values of the 
Casimir force within some definite separation region, almost unchanged. The stability of 
the Casimir force opens new opportunities to use the Casimir effect as a test for long-range 
hypothetical interactions and for the diagnostic purposes. For example, some kind of the 
inverse problem could be utilized, i.e. the measured force- distance relations be exploited to 
determine the fundamental characteristics of solids (such as the plasma frequency). 
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FIG. 1: 15 x 15 /um 2 atomic force microscope image of the Au coating on the plate. The topography 
of the coating on the sphere is similar. 
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FIG. 2: Typical cross sections of the atomic force microscope image of the Au coating on the plate 
with (a) constant x and (b) constant y. 



28 



Tables 



29 





TABLE I: Fractions Vi of the surface area covered by rouj 


diness with heights hi. 


i 


/i,- (nm) 


Vi 


1 





1.06 x 10~ 3 


2 


1 


5.086 x 10~ 2 


3 


2 


0.33511 


4 


3 


0.45863 


5 


4 


0.13695 


6 


5 


1.586 x 10~ 2 


7 


6 


1.24 x 10~ 3 


8 


7 


1.6 x 10~ 4 


9 


8 


4 x 10~ 5 


10 


9 


2 x 10~ 5 


11 


10 


1 x 10" 5 


12 


11 


1 x 10~ 5 


13 


12 


1 x 10~ 5 


14 


13 


1 x 10~ 5 


15 


14 


1 x 10~ 5 


16 


15 


1.2 x 10~ 5 


17 


16 


8 x 10~ 6 
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TABLE II: Corrections factors to the ideal Casimir force at different separations due to finite con- 
ductivity (r) c ), surface roughness (rj r ) and both finite conductivity and surface roughness (r] C;r and 
r]™ r in the method of the geometrical averaging and in the multiplicative approach, respectively). 

z = 62 nm z = 70 nm z = 80 nm z = 90 nm 

0.4681 0.4964 0.5218 

1.0017 1.0013 1.0010 

0.4687 0.4669 0.5223 

0.4689 0.4670 0.5223 



Vr 



m 
'Ic.r 



0.4430 
1.0022 
0.4436 
0.4440 



31 



